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Abstract 

In studies of the dynamics of charged particles in a cyclotron magnetic field, the specified field is 
generally B z in the z — midplane where B r and Bg are zero. B r (r,9,z) and Bg(r,8,z) are simple to 
determine through a linear expansion which assumes that B z is independent of z. But, an expansion to 
only first order may not be sufficient for orbit simulations at large z. This paper reviews the expansion 
of a specified B z (r,8, z — 0) out of the z — midplane to arbitrary order, and shows simple examples 
worked out to 4th order. 
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1 Introduction 

Many programs and software packages simulate orbits of charged particles in magnetic fields as noted in 
Table 1. Some use a single reference orbit and represent the magnetic field as a series of matrices, one for 
each magnetic element. Other programs allow more realistic spiral orbits by using magnetic field equations 
or interpolated magnetic field maps for point to point Runge-Kutta track integrators. Some programs such 
as COSY Infinity use exact analytical determination of magnetic field derivatives; others use numerical dif- 
ferentiation. An important preliminary step before simulation of fields produced by magnets is investigation 
of orbit dynamics in an idealized magnetic field which satisfies Maxwell's equations in vacuum to sufficient 
order. Magnetic fields for cyclotrons [T| are most conveniently expressed in cylindrical (r,6,z) coordinates in 
which the z direction is the cyclotron rotation axis. In many applications, an ideal field is specfied in the 
cyclotron median plane which is the z — plane where B r and Bg vanish. Determination of (B r ,Bg,B z ) 
for z 7^ to sufficient order (beyond 1st order) is necessary for detailed orbit and dynamics studies [3J. 
Reduction of the phase space of a muon beam through ionization cooling in an inverse cyclotron is being 
researched [3] as part of the muon R&D for possible future facilities such as a neutrino factory and muon 
collider g]. 

2 Expansion of Field Out of the Midplane 

The importance of correct field expansion for reliable particle tracking simulation and beam dynamics analysis 
has been recognized, and the expansion around the reference orbit of a bent-solenoid geometry has been 
calculated [14] . Expansions to eighth order are done with the aid of Mathematica. A static magnetic field 
can be expanded out of a plane because it must satisfy the source- free Maxwell equations in vacuum, V • B = 
and V x B = 0. Only dB x /dy — dB y /dx constrains fields in the plane, so any function can be used for B z . 
One gets the three field components out of the plane as follows 

dB z _ dB x dB y dB x _ dB z dB y _ dB z 

dz dx dy dz dx dz dy 

The expansion of a magnetic field specified in a plane has been addressed by Zgoubi for cartesian coor- 
dinates [TS]. Zgoubi uses Taylor expansion with V • B = and V x B = 0. Here we work out a solution in 
cylindrical coordinates to arbitrary order. Median plane antisymmetry is assumed so that 
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Table 1: Computer programs for tracking charged particles in magnetic fields and designing accelerators. 
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A rigorous proof of these expansions is provided in the final section of this paper. 

In most implementations of simulating a magnetic field, the expansions of (B r , Bg, B z ) cannot be ex- 
pressed in closed form. In such cases, only a finite number of terms of the magnetic field expansion can be 
used. The order of the expansion is defined as the highest power of z in the summation. To odd (even) order 
to, B r and Bg have ( y) terms while B z has terms. When B is expanded to odd order, 

VxB = and V • B ^ 0; when B is expanded to even order, V ■ B = and V x B / 0. An even order 
expansion with V • B = is needed to satisfy the Hamiltonian and allow a symplectic map |16j . 
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3 Example Fields to 4th Order 



This section shows some example midplane fields B z (r, 9, z = 0) expanded to 4th order to obtain B r (r, 9, z), 
Bg(r,9 7 z) 7 and B z (r,6,z). These example fields have fairly simple radial and azimuthal dependences and 
were checked with Mathemtica. The 4th order expansion of a more complex field, one with a constant radial 
field index, k and spiral angle, a, such as B z (z = 0) = cr fe (l — f sin[N(8 — tanaln(r/ro))]), yields pages of 
output so that the use of a program such as Mathematica is necessary. Mathematica includes a command 
which produces output in FORTRAN or C syntax which can be copied directly into a routine which generates 
a grid of magnetic field points. Example Mathematica and FORTRAN files are presented in Appendices A and 
B. 

3.1 B Only Dependent on Radius (No Sectors) 

Is this case, B(r, 6, z = 0) is B(r, z = 0) so that 

dB z 3 d fd 2 B 1 dB 

B r — Z— — — — 



dr 6 dr \ dr 2 r dr 



B e = 



. z 2 fd 2 B ldB\ z 4 ( d 4 B d 2 fldB\ ld 3 B Id fldB 
B z = B{r, z = 0) 



2 \ dr 2 r dr J 24 \ dr 4 dr 2 \r dr J r dr 3 r dr \r dr 

An example of such a field is a non-sectored field in the middle of Hclmoltz coils. At z = 0, the B z is 
maximal at r — and decreases very slowly with radius so that dB/dr and higher derivatives are small, but 
still large enough to provide weak focusing in r and z. 

3.2 Azimuthal Field Sectors: No Spiral Angle nor Midplane Radial Dependence 

For this field B(r, 9,z = 0) is B(9, z = 0): an example is B = B (l - fsm(N6)) so that 

B r = l(^Y (fB N 2 sm(N9)) 

Be = -^(fB Ncos(N9)) - 1 ^ (fB N 3 cos(N9)) 

B z = B (l - fsm(N9)) - l - ^ (fB N 2 sin(N9)) + ^ (jj) * (fB N 2 sm(N9))(4 - N 2 ) 
A sectored field such as this provides stronger focusing in z than a non-sectored field. 

3.3 Azimuthal Field Sectors: Constant Spiral Angle, No Radial Field Index 

For this field, B(r, 9, z = 0) is B(9, z = 0): a field with constant spiral angle a is 

B = B Q (l - fsm[N(6 - tanaln(r/r ))]) so that 

B r = - (fB N tan a) (cos[A (9 - tanaln(r/r ))])+ 
r 

\ (-Y (fB N 2 sec 2 a)(Nt&nacos{N(9 - tanalog(r/r ))] + 2sin[iV(6l - tan a log(r/r ))]) 
6 \rJ 
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B z = B (l - fsm[N(6 - tana ln(r/r ))])- 
\ (-) 2 (/^o^ 2 sec 2 a)(sinLV(0 - tan a ln(r/r ))])+ 

^ (^) 4 (f B ° N2 se ° 2 a ) { ~ N2 se ° 2 a $HN(0 - tan a ln(r/r ))] + 
4(sin[iV(fl - tanaln(r/ro))] + 7Vtanacos[iV(6> - tanaln(r/r ))]} 

The spiraled magnetic field provides stronger focusing in z than the non-spiraled sectored field. When a 
= 0, the magnetic field components reduce to those of the previous azimuthally sectored field with no spiral 
angle nor midplane radial dependence. 

3.4 Azimuthal Field Sectors: No Spiral Angle, Constant Radial Field Index 

r dB z (z = 0) 
~ B z (z = 0) dr 

For this field, B(r,9,z = 0) = cr fc (l - fsm(N6)). Then 

B r = (^) (cr fc )fc(l - fsm(N6)) - ^ (^f (cr k )(k - 2)(fN 2 sm(N6) + fc 2 (l - fsm(N6))) 

B g = - [^j (cr k )fNcos(N0) - ^ ^ (cr k )(fN cos(N9))(N 2 - k 2 ) 
B z = cr k (l- fsm(N9))- 



2 (~) (cr k )(fN 2 sm(N9) + k 2 (l - f sm(N9)))+ 
(J?)* \cr k ) [fN 2 sin(N9)(2(k 2 - 2k + 2) - iV 2 ) + (k - 2) 2 fc 2 (l - / sin(JV0))] 



1 p\ 4 
24 Vr 



4 Proof of Expansion Out of the Plane 

Let 

B(r, 9, z) = fB r {r, 9, z) + 9B g (r, 9, z) + zB z (r, 9, z) = B ± (r, 9, z) + zB z (r, 9, z), 

B ± (r,9,z) = j:Z A k (r,9)z k , 
B z (r,9,z) = J2Z B k (r,9)z k 1 

For any vector field V and scalar field S, 

It is necessary and sufficient to show that the specification of 
B (r, 0), V • B = 0,VxB = 04 
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(Vx • A fc ) + (fc + l)B k+l = => 

S fe +i = -^(Vi • A fe ) for fc = 0, 1, 2, . . . 

V x B = => 

(v±+z^ x (Bx + = 
(Vx x B ± ) + (Vx x (zB z )) + x B ± ) = 

(Vx x Bx) - z x (V±B Z ) + z x (2fL) = 



(Vx x Bx) - 5 x (VxB z - 



0. 



(Vx x Bx) is parallel to z and £ x |^Vii? z — ^§^-^j is perpendicular to z. 



Therefore, (Vx x Bx) = and (y±B z - ?§±) = 0. 
Then (VxS, - %-) = => EZo(V±B k )z k - £~ i fcAfcZ^ 1 = 
J2Zo ((VxSfe) - (fc + l)A fe+1 ) z k = => 
Afc+i = ^(ViBfc) for fc = 0, 1, 2, . . . 



Bk+i = -^(Vx • A fc ) and A fe+1 = ^(V \B k 
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fc = ^ 
Ai = V±B 
B2 = —\V\Bq 

A = 
Bx = 

fc = 1 

A 2 = ±V ± Bi = 

B 3 = -(37(2) V i Bl = 



A 3 = iVxS 2 = -(3J(2) Vx (Viflo) = -ijVx (viBo) 

B 4 = -(4jp) vis 2 = (4)( 3 )(2) vi (v^So) = irVi (vis ) 



Continuing recursively, 

fc = 2n - 1 

A 2n = for n > 1 
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£?2n-i = for n > 1 
k = 2n => 

A 2«+i = S&V ± ((V 2 )^ ) for n > 1 
fi 2n = ^r(V 2 r|5 for n > 1 

The proof is complete. 

5 Summary 

The procedure for expanding a specified B z (r, 9, z — 0) out of the z = midplane is described. Examples of 
expansions out to 4th order are included. Fig. 1 shows an example of a magnetic field. 



Bz(x,y,z=0) [T] 




Figure 1: B z (x,y,z = 0) = 1.77r a6 (l + V^cos (60)) (1 + tanh (25(r - 0.2))/2 including the outer 6-sector 
focus field which merges into the inner magnetic bottle. B z (x,y,z = 0) ranges from -0.6 to 3.2 Tesla. The 
inner bottle field is generated by two r = 0.2 m circular coils, located 0.2 m above and below the midplane. 
The B z field at r = 0, z = is 2.4 Tesla. 
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ren, M. Berz, A. Bogacz, and R. Palmer for their advice and help. 
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A Appendix: Mathematica Example 

An example of a Mathematica file which evaluates the coefficients of B z (r,9,z = 0) = cr k (l + f cos(Nd)) 
up to 4th order is shown. This B z (r,9,z = 0) is simpler than the field of Fig. Q] so that the example file 
and output are reasonably short and clear. The example file also formats the expressions of the coefficients 
to FORTRAN. The next appendix is an example FORTRAN routine which uses the output FORTRAN expressions 
from this file to generate a grid of magnetic field points used by a G4Bcamline input file (Appendix C). This 
Mathematica file shows the output of each Mathematica command indented. 

(* Set b as Bz(z = 0) with (r, t) as radius and theta. *) 
b[r_,t_] = (c r~k)(l+f Cos[n ( t)]) 

c r~k (1+f Cos [n t] ) 

(* Set trans as d~2b/dr~2 + (l/r)db/r + (l/r~2)d~2b/dt~2 . *) 

trans[r_,t_] = D[b[r,t] ,{r,2>] + (l/r)D [b [r,t] , {r, 1>] + (l/r~2)D[b[r,t] ,{t,2>] 

-c f n~2 r~(-2+k) Cos[n t]+c k r~(-2+k) (1+f Cos [n t])+c (-1+k) k 
r~(-2+k) (1+f Cos[n t] ) 

(* FullSimplif y does an algebraic simplification *) 
FullSimplifyra 

c r~(-2+k) (k~2+f (k-n) (k+n) Cos [n t] ) 

(* trans2 is a double application of the derivatives of trans. *) 
trans2 [r_,t_] = D [trans [r ,t] , {r ,2}] + (l/r)D [trans [r,t] , {r, 1}] + 
(l/r~2)D [trans [r,t] ,{t,2>] 

-c f (-3+k) (-2+k) n~2 r~(-4+k) Cos [n t]+c (-3+k) (-2+k) k r~(-4+k) 
(1+f Cos[n t])+c (-3+k) (-2+k) (-1+k) k r~(-4+k) (1+f Cos [n t])+ 



<S 



(-c f k n~2 r~(-2+k) Cos [n t] -c f (-1+k) k n~2 r~(-2+k) Cos [n t] + 

c f iT4 r~(-2+k) Cos [n t])/r'2+l/r (-c f (-2+k) n~2 r~(-3+k) Cos [n t] + 

c (-2+k) k r~(-3+k) (1+f Cos [n t])+c (-2+k) (-1+k) k r~(-3+k) 

(1+f Cos[n t])) 

FullSimplifym 

c r~(-4+k) ((-2+k) ~2 k~2+f (-2+k-n) (k-n) (-2+k+n) (k+n) Cos[n t] ) 

(* Set field components, apply simplifications, convert to *) 
(* FORTRAN format. *) 

(* Use ol, o2, o3, and o4 as switches to set order of expansion. *) 
(* For 2nd order, set ol = 1, o2 = 1 , o3 = o4 = 0. *) 

(* Set the 0th, 2nd, and 4th order terms of Bz(z). *) 
bz0[r_,t_] = b[r,t] 

c r~k (1+f Cos[n t] ) 
FullSimplifyra 

c r~k (1+f Cos[n t] ) 
FortranForm [°/ ] 

c*r**k*(l + f*Cos(n*t)) 

bz2[r_,t_] = -(o2) (z~2/2)trans[r,t] 

-(1/2) o2 z~2 (-c f n~2 r"(-2+k) Cos[n t]+c k r~(-2+k) (1+f Cos [n t])+ 
c (-1+k) k r~(-2+k) (1+f Cos [n t])) 

FullSimplify[7 ( ] 

-(1/2) c o2 r~(-2+k) z~2 (k~2+f (k-n) (k+n) Cos[n t] ) 

FortranForm [°/ ] 

-(c*o2*r**(-2 + k)*z**2*(k**2 + f*(k - n)*(k + n) *Cos (n*t) ) ) /2 . 

bz4[r_,t_] = (o4) (z~4/24)trans2 [r,t] 

1/24 o4 z~4 (-c f (-3+k) (-2+k) n~2 r~(-4+k) Cos [n t]+ 

c (-3+k) (-2+k) k r~(-4+k) (1+f Cos [n t])+ 

c (-3+k) (-2+k) (-1+k) k r~(-4+k) (1+f Cos[n t])+ 

(-c f k n~2 r~(-2+k) Cos [n t]-c f (-1+k) k n~2 r"(-2+k) Cos [n t]+ 

c f n~4 r~(-2+k) Cos [n t])/r"2+l/r (-c f (-2+k) n"2 r"(-3+k) Cos [n t]+ 

c (-2+k) k r~(-3+k) (1+f Cos [n t])+ 

c (-2+k) (-1+k) k r~(-3+k) (1+f Cos [n t]))) 

FullSimplify[7 ( ] 

1/24 c o4 r"(-4+k) z"4 ((-2+k) ~2 k"2+f (-2+k-n) (k-n) (-2+k+n) (k+n) 
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Cos [n t] ) 

FortranForm [%] 

(c*o4*r**(-4 + k)*z**4* 
((-2 + k)**2*k**2 + 

f*(-2 + k - n)*(k - n)*(-2 + k + n)*(k + n) *Cos (n*t) ) ) /24 . 

(* Set the 1st and 3rd order terms of Br(z) . *) 
brl[r,t] = (ol)z D[b[r,t],{r,l>] 

c k ol r~(-l+k) z (1+f Cos[n t] ) 
FullSimplifyra 

c k ol r~(-l+k) z (1+f Cos[n t] ) 
FortranForm [°/«] 

c*k*ol*r**(-l + k)*z*(l + f*Cos(n*t)) 

br3[r,t] = -(o3) (z~3/6)D [trans [r ,t] , {r, 1}] 

-(1/6) o3 z~3 (-c f (-2+k) n~2 r~(-3+k) Cos[n t]+ 

c (-2+k) k r~(-3+k) (1+f Cos [n t])+ 

c (-2+k) (-1+k) k r~(-3+k) (1+f Cos [n t])) 

FullSimplify[7,] 

-(1/6) c (-2+k) o3 r~(-3+k) z~3 (k~2+f (k-n) (k+n) Cos [n t] ) 

FortranForm [%] 

-(c*(-2 + k)*o3*r**(-3 + k)*z**3* 

(k**2 + f*(k - n)*(k + n) *Cos (n*t) ) ) /6 . 

(* Set the 1st and 3rd order terms of Btheta(z) . *) 
btl[r,t] = (ol)z (1/r) D[b[r,t] ,{t,l>] 

-c f n ol r~(-l+k) z Sin[n t] 
FullSimplifyra 

-c f n ol r~(-l+k) z Sin[n t] 
FortranForm [%] 

-(c*f *n*ol*r**(-l + k)*z*Sin(n*t)) 

bt3[r,t] = -(o3) (z~3/6) (l/r)D[trans[r,t] ,{t,l>] 

-1/(6 r) o3 z~3 (-c f k n r~(-2+k) Sin[n t] - 

c f (-1+k) k n r~(-2+k) Sin[n t]+c f n~3 r~(-2+k) Sin[n t] ) 
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FullSimplifyra 

1/6 c f (k-n) n (k+n) o3 r~(-3+k) z~3 Sin[n t] 
FortranForm [°/,] 

(c*f*(k - n)*n*(k + n)*o3*r**(-3 + k)*z**3*Sin(n*t))/6. 



B Appendix: FORTRAN Example 

This program takes the output expressions of B z {z), B r (z), and Bq(z) from Mathcmatica and prints out a 
grid of magnetic field values recast as B x , B y , and B z . The output file generated by this routine is fort. 71. 
With further formatting described in the initial comments of the routine, the output file can be used by a 
G4Bcamline input file to generate a magnetic field for an anticyclotron simulation. In this example, the field 
componenets are expanded to 4th order in z. 

PROGRAM FOURTH_ORDER_ARXIV_EXAMPLE 

C 

C Prints out grid of (Bx, By, Bz) in (x, y, z) for a file suitable for the 

C G4Beamline command 'f ieldmap' . 

C 

C The output of this program is fort. 71. 
C 

C Distances in the fort. 71 file are in mm., and magnetic fields are 

C in tesla. (Distances in this FORTRAN routine are in meters.) 

C 

C The following must be added to fort. 71 for the file to work in G4Beamline: 
C 

C grid X0= [smallest x] Y0= [smallest y] Z0= [smallest z] nX= [number of x points] 
C nY= [number of y points] nZ= [number of z points] dX=[dist. between x points] 
C dY=[dist. between y points] dZ=[dist. between z points] 
C data 
C 

IMPLICIT NONE 
INTEGER I.J.L 
REAL*8 x,y,z,r,t 

REAL*8 bzO , bz2 , bz4 , br 1 , br 3 , bt 1 , bt 3 
REAL*8 bx,by,bz,br ,bt 

REAL*8 PI,E 
REAL*8 b,c,k,s,f ,n 

REAL*8 ol,o2,o3,o4 

E = 2.718281828d+0 
PI = AC0S(-1.0d+0) 

C Bz(z=0) = (c*r~k)*(l+f *cos(n(t))) 
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Bz(z = 0) parameters 
b = 25.0d+0 
c = 1.77d+0 
k = 0.6d+0 

s = (PI/180. 0d+0)*0.0d+0 
f = sqrt(2.0d+0) 
n = 6.0d+0 



C ol through o4 set the order of the expansion in z of Bz(z=0) 



c 


4th 


order : 


o4 = 


o3 


= o2 


= ol 


= 1 




c 


3rd 


order : 


o4 = 


0, 


o3 = 


o2 = 


ol = 


1 


c 


2nd 


order : 


o4 = 


o3 


= o, 


o2 = 


ol = 


1 


c 


1st 


order : 


o4 = 


o3 


= o2 


= o, 


ol = 


1 


c 


0th 


order : 


o4 = 


o3 


= o2 


= ol 


= 





ol = 


1 


,0d+0 


o2 = 


1 


.0d+0 


o3 = 


1 


.0d+0 


o4 = 


1 


,0d+0 



C Loop over x, y, z 

DO 100 I = 1, 201 
DO 200 J = 1, 201 
DO 300 L = 1, 61 

C x, y, z, r in meters 

C t in radians 

C br, bt, bz in Tesla 



x = -1.01d+0 + 0.01d+0*I 
y = -1.01d+0 + 0.01d+0*J 
z = -0.31d+0 + 0.01d+0*L 

r = sqrt(x*x+y*y) 

t = atan2(y,x) 



C br, bt, bz formula determined through and copied and pasted from Mathematica 



C bz terms up to 4th order 

bzO = c*r**k*(l + f*Cos(n*t)) 



bz2 = - (c*o2*r** (-2 + k) *z**2* (k**2 + 
f*(k - n)*(k + n)*Cos(n*t)))/2. 



bz4 = (c*o4*r**(-4 + k)*z**4* 
((-2 + k)**2*k**2 + 

f*(-2 + k - n)*(k - n)*(-2 + k + n)*(k + n) *Cos (n*t) ) ) /24 . 



C br terms up to 4th order 

brl = c*k*ol*r**(-l + k)*z*(l + f*Cos(n*t)) 
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br3 = -(c*(-2 + k)*o3*r**(-3 + k)*z**3* 
(k**2 + f*(k - n)*(k + n) *Cos (n*t) ) ) /6 . 

C bt terms up to 4th order 

btl = -(c*f * n *ol*r**(-l + k) *z*Sin(n*t) ) 

bt3 = (c*f*(k - n)*n*(k + n)*o3*r**(-3 + k) *z**3*Sin(n*t) ) /6 . 

C Add each order to get bz, br, bt 
bz = bzO + bz2 + bz4 

br = brl + br3 

bt = btl + bt3 

C Transform br, bt to bx, by 

bx=br*C0S (t) -bt*SIN (t) 

by=br*SIN(t)+bt*COS(t) 

WRITE(71,174) 1000*x,1000*y,1000*z,bx,by,bz,0.0 

174 FORMAT (F8 . 1 , " , " ,F8 . 1 , " , " ,F8 . 1 , " , " ,F9 . 3 , " , " ,F9 . 3 , " , " ,F9 . 3 , 

1 ", 0.0, 0.0, \F7.3) 

300 CONTINUE 
200 CONTINUE 
100 CONTINUE 

END 

Here are the first 10 lines in the fort. 71 output file. 



head fort. 71 



-1000. 


.0, 


-1000. 


o, 


-300. 


0, 


-3. 


,318, 


3. 


,712, 


2. 


162, 


0. 


.0, 


0. 


o, 


0. 


.000 


-1000. 


.0, 


-1000. 


.0, 


-290. 


.0, 


-3. 


.160, 


3. 


,542, 


2. 


.163, 


0. 


.0, 


0. 


.0, 


0. 


.000 


-1000 


,0, 


-1000, 


o, 


-280. 


0, 


-3. 


.007, 


3. 


.375, 


2. 


.164, 


0. 


.0, 


0. 


.0, 


0. 


,000 


-1000 


.0, 


-1000. 


.0, 


-270. 


.0, 


-2 


.859, 


3. 


,214, 


2. 


.165, 


0. 


.0, 


0. 


.0, 


0. 


.000 


-1000 


,0, 


-1000. 


o, 


-260. 


0, 


-2. 


.716, 


3. 


,057, 


2. 


.166, 


0. 


.0, 


0. 


.0, 


0. 


.000 


-1000 


.0, 


-1000. 


.0, 


-250. 


,0, 


-2 


,576, 


2 


,904, 


2. 


.167, 


0. 


.0, 


0. 


o, 


0. 


.000 


-1000. 


• 0, 


-1000. 


,0, 


-240. 


0, 


-2 


.441, 


2. 


.756, 


2. 


.168, 


0. 


o, 


0. 


o, 


0. 


.000 


-1000. 


.0, 


-1000. 


.0, 


-230. 


.0, 


-2. 


.310, 


2. 


.611, 


2. 


.169, 


0. 


.0, 


0. 


o, 


0. 


.000 


-1000. 


• 0, 


-1000. 


o, 


-220. 


0, 


-2. 


.182, 


2. 


.471, 


2, 


170, 


0. 


o, 


0. 


o, 


0. 


.000 


-1000. 


.0, 


-1000. 


o, 


-210. 


.0, 


-2. 


.058, 


2. 


.334, 


2. 


.171, 


0. 


.0, 


0. 


o, 


0. 


.000 



With a grid spacing of 10 mm ranging from —1000 mm < (x,y) < 1000 mm and —300 mm < z < 
300 mm, the following lines need to be placed at the start of fort .71 to match the syntax required by the 
G4Beamline command f ieldmap 

param normB=1.0000 normE=0.0000 

grid X0=-1000 Y0=-1000 Z0=-300 nX=201 nY=201 nZ=61 dX=10 dY=10 dZ=10 
data 
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C Appendix: G4Beamline Example 



An example G4Beamline [5] input file of positive muons being decelerated in an anticycloton is shown. The 
file test_bfield_6sector_rk06_0d_2nd_order_cos_example.dat as specified in the fieldmap field_grid 
f ile=test_bf ield_6sector jrk06_0d_2nd_order_cos_example . dat command is needed to generate the outer 
sectored focusing field. The inner magnetic bottle is generated by the field from two current carrying coils. 
Muons lose energy by passing through 1.0° lithium hydride wedges. The uniform 1.0° wedge angle in this 
example file is simpler than the moderator configuration used for the anticyclotron. The anticyclotron 
includes an inner cylinder of 0.1 bar helium and six LiH wedges whose thickness decreases adiabatically from 
4.2 mm at r = 550 mm to 0.008 mm at r = 65 mm. The output is a root file, AllTracks . root, which 
shows all the kinematic information of each muon along its orbit, and which can be opened by Historoot. 

# example G4Beamline file for anticyclotron simulation, April 11, 2011 
* 

# use TJR recommended physics routines 
physics QGSP_BERT doStochastics=l disable=Decay 

# reference particle with no stochastic processes 
particlecolor ref erence=l , 1 , 1 

reference ref erenceMomentum=180 particle=mu+ \ 

beamX=0.0 beamY=500.0 beamZ=0.0 meanXp=0.0 meanYp=0.0 rotation=X90 , Z90 

# simulate 3 muons with stochastic processes (since doStochastics=l) 
beam gaussian particle=mu+ nEvents=3 \ 

beamX=0.0 beamY=500.0 beamZ=0.0 meanXp=0.0 meanYp=0.0 \ 
sigmaX=0 sigmaY=0 sigmaXp=0 sigmaYp=0 \ 

meanMomentum=180 . 00 sigmaP=0 meanT=0 sigmaT=0 rotation=X90 , Z90 



# keep only positive muons 
trackcuts keep=mu+ maxTime=200000 . 



# output file for Historoot showing kinematic information for orbits 

trace nTrace=3 format=root f ilename=AHTracks . root oneNTuple=l primary0nly=l 

# output file for Historoot showing kinematic information when particles are lost 
beamlossntuple loss_ntuple format=ascii 



param maxStep=5. SteppingVerbose=l 

# lithium hydride, density (0.82 g/cnT3) from Wikipedia 
material lih z=4 a=8 density=0.82 state=s 

# Define LiH wedges, each with 1 degree angle 

tubs gascylinder_l innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=30-0.5 f inalPhi=30+0 . 5 color=l,0,0 material=lih 
tubs gascylinder_2 innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=90-0.5 f inalPhi=90+0 . 5 color=l,0,0 material=lih 
tubs gascylinder_3 innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=150-0.5 f inalPhi=150+0 . 5 color=l,0,0 material=lih 
tubs gascylinder_4 innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=210-0.5 f inalPhi=210+0 . 5 color=l,0,0 material=lih 
tubs gascylinder_5 innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=270-0.5 f inalPhi=270+0 . 5 color=l,0,0 material=lih 
tubs gascylinder_6 innerRadius=2 outerRadius=600 length=200 \ 

initialPhi=330-0.5 f inalPhi=330+0 . 5 color=l,0,0 material=lih 
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# Define coils of inner magnetic bottle 
coil coill \ 

innerRadius=195 \ 
outerRadius=205 \ 
length=50 \ 
material=Vacuum \ 
f ilename=coill . dat 

coil coil2 \ 
innerRadius=195 \ 
outerRadius=205 \ 
length=50 \ 
material=Vacuum \ 
f ilename=coil2 . dat 

# Set currents for magnetic bottle coils which provide inner bottle magnetic field 
solenoid cl \ 

coilName=coill \ 
current=(0.2)*10803.6 \ 
color=0.3,l,0 \ 
alternate=0 

solenoid c2 \ 

coilName=coill \ 
current=(0.2)*10803.6 \ 
color=0.3,l,0 \ 
alternate=0 

# Grid of magnetic field points for outer sectored field 

f ieldmap f ield_grid f ile=test_bf ield_no_poly_bottle_6sector_rk06_0d_2nd_order_beq25_cos . dat 

# Place wedges, coils, and magnetic field map grid 
place gascylinder_l z=0 

place gascylinder_2 z=0 
place gascylinder_3 z=0 
place gascylinder_4 z=0 
place gascylinder_5 z=0 
place gascylinder_6 z=0 

place cl x=0 y=0 z=200 
place c2 x=0 y=0 z=-200 

place field_grid x=0 y=0 z=0 

Fig.[2]shows a Historoot visualation of the orbits of three muons in this G4Beamline example simulation. 
The muons stop in about 125 ns. 

D Appendix: ICOOL Example 

Another program which can simulate an anticyclotron with lithium hydride wedges is ICOOL [7]. The 
coordinate system of ICOOL is an (x, y, z) coordinate system in which the z coordinate is along the beamline. 
This can be converted to an (r, 9, z) system by incorporating a bent solenoid so that x becomes — r, y becomes 
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Figure 2: Orbits of three color coded unions in an G4Beamline simulated anticyclotron. This output is 
generated by Historoot. The orbits show the transition between the outer sectored field and the inner 
azimuthally symmetric magnetic bottle. 



z, and the reference path coordinate, z, becomes 0. ICOOL does not have the visualization capabilites of 
G4Beamline so software to check element placement, orbits, and dynamics must be added. ICOOL can do 
its own magnetic field expansion [14 to check an orbit at a particular radius. This can provide a cross check 
of G4Beamline. 

ICOOL can be run on UNIX, PC, and Macintosh platforms. The command file needs to be called 
for001.dat. In this particular application of ICOOL, two more input files are used: for003.dat which 
defines each particle of the beam and for056.dat which defines the magnetic field. The output file in the 
example simulation is for009.dat which lists the kinematic information of each particle. Postprocessing of 
this file is needed to visualize the orbits and to evaluate the dynamics. 

Information about ICOOL including a User's Guide and Reference Manual can be found at 

https : //pubweb .bnl . gov/~f ernow/ icool/v326/ 

The next sections show example input files which generate an anticyclotron simulation when the executable 
icool from 

https : //pubweb .bnl . gov/ ~f ernow/ icool/v326/icool . exe 
is run. 



D.l input for001.dat command file 

This command files uses a bent solenoid magnetic field to simulate a cyclotron field. The BSDL command 
sets the initial muon momentum to 0.18 GeV/c, the inverse radius of the bent solenoid to 2.0 m _1 , and reads 
the amplitude, period, and initial offset of the sinusoidal vertical magnetic field (B y (z) where y is vertical 
and z is along the reference path) from the file fort. 56. In this ICOOL example, muons lose energy by 
passing through 1.0° lithium hydride wedges. The magnetic field in this example does not vary with radius 
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is in the G4Beamline example of Appendix C. The geometry of the anticyclotron is specified by sregions. 
One sregion of the lithium hydride wedge is followed by 10 sregions of vacuum, and the sequence of 11 
sregions is repeated 600 times. The total length of the 11 sregions is 0.523598775 m which corresponds 
to 1/6 of the circumference of a 0.5 to circle. This arrangement is set for a 6-scctor anticyclotron with a 
bent solenoid radius of 0.5 to. 

test ring 180 MeV (H5b) 

$cont npart=200 nsections=l varstep= . false . nprnt=-3 prlevel=l 
ntuple= . false . rtuple= . false . rtuplen=l output 1= . true . phasemodel=3 
fsav=. false, f savset= . true . izfile=140 bgen=. false. $ 

$bmt $ 

12 1.1 !ntyp typ(le 2mu 3pi 4K 5p) frac dist(l=gauss 2=ring) 

-0.010 0.019 0. 0. 0. 0.203 ! means: x y z px py pz 

0.0425 0.0425 .18 .030 .030 0.024 Isigmas: 

1 ! correlations 

2 .19 .40 ! palmer GeV/A~2 beta 

$ints ldecay=. false. declev=l ldedx=.true. lstrag= . true . lscatter= . true . 
delev=2 straglev=4 scatlev=4 $ 

$nhs $ 
$nsc $ 
$nzh $ 

$nrh nrhist=0 $ 
$nem $ 
$ncv $ 

SECTION 

REFP 

2 .18 0. 0. 3 !typ refp tO gradO mode (3=const p 4=with acc) 

DENS 
LIH 1 
BEGS 



CELL ! regular cell 

1 

. FALSE . 

BS0L ! multipole field input 

4. 56 .180 51111111111 2.0 ! mode, file, momentum, order, 1/r switch, 10 scale-factors, 1/: 



REPEAT 
1200 

OUTPUT 

SREGION ! define a region 
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.023598775 1 0.0001 ! length, 1 radial subregion, st< 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

LIH VAC ! LiH material surrounded by vacuum 

WEDGE 

1 0.5 0.01 180 1.0 1.5 0. 0. 0. 0. ! angle RVERT ZVERT Azymuth Apx 

OUTPUT 

SREGION ! define a region 

.05 1 0.001 ! length, 1 radial subregion, step 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

VAC ! vacuum 

CBLOCK ! cylindrical block geometry 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

OUTPUT 

SREGION ! define a region 

.05 1 0.001 ! length, 1 radial subregion, step 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

VAC ! vacuum 

CBLOCK ! cylindrical block geometry 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

OUTPUT 

SREGION ! define a region 

.05 1 0.001 ! length, 1 radial subregion, step 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

VAC ! vacuum 

CBLOCK ! cylindrical block geometry 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

OUTPUT 

SREGION ! define a region 

.05 1 0.001 ! length, 1 radial subregion, step 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

VAC ! vacuum 

CBLOCK ! cylindrical block geometry 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

OUTPUT 

SREGION ! define a region 

.05 1 0.001 ! length, 1 radial subregion, step 

1 0. 2 ! radial extent 

NONE ! no associated field 

0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 
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VAC 
CBLOCK 
0. 0. 0. 



0. 



! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 0. 



OUTPUT 
SREGION 
.05 1 0.001 
1 0. 2 
NONE 

0. 0. 0. 0. 
VAC 

CBLOCK 

0. 0. 0. 0. 

OUTPUT 
SREGION 
.05 1 0.001 
1 0. 2 
NONE 

0. 0. 0. 0. 
VAC 

CBLOCK 

0. 0. 0. 0. 

OUTPUT 
SREGION 
.05 1 0.001 
1 0. 2 
NONE 

0. 0. 0. 0. 
VAC 

CBLOCK 

0. 0. 0. 0. 

OUTPUT 
SREGION 
.05 1 0.001 
1 0. 2 
NONE 

0. 0. 0. 0. 
VAC 

CBLOCK 

0. 0. 0. 0. 

OUTPUT 
SREGION 
.05 1 0.001 
1 0. 2 
NONE 

0. 0. 0. 0. 
VAC 

CBLOCK 

0. 0. 0. 0. 



! define a region 
! length, 1 radial subregion, step 
! radial extent 
! no associated field 
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 



! define a region 
! length, 1 radial subregion, step 
! radial extent 
! no associated field 
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 



! define a region 
! length, 1 radial subregion, step 
! radial extent 
! no associated field 
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 



! define a region 
! length, 1 radial subregion, step 
! radial extent 
! no associated field 
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 



! define a region 
! length, 1 radial subregion, step 
! radial extent 
! no associated field 
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! vacuum 

! cylindrical block geometry 

0. 0. 0. 0. 0. 
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ENDREPEAT 

ENDCELL 
ENDSECTION 



D.2 input for003.dat beam file 

This input file defines the initial positions and momenta of three input muons. The initial positions are set 
to zero with respect to the 0.5 m reference radius of the bent solenoid field. The initial momentum of the 
three muons is 0.18 GeV/c in the z direction. 



title 

00000000 

10 2 0.0000E+00 1.0 

0.00E-00 0.0000E+00 0.0000E+00 0.0000E-03 0.0000E+00 1.8000E-01 

2 2 0.0000E+00 1.0 

0.00E-00 0.0000E+00 0.0000E+00 0.0000E-03 0.0000E+00 1.8000E-01 

3 2 0.0000E+00 1.0 

0.00E-00 0.0000E+00 0.0000E+00 0.0000E-03 0.0000E+00 1.8000E-01 



D.3 input for056.dat field file 

This file sets the parameters of the bent solenoid field. This field is sinusoidal with respect to the z coordinate 
(In ICOOL, the vertical component is y and z is the coordinate along the reference path.). The period is 
about 0.5236 m corresponding to 1/6 the circumference of a circle with reference radius of 0.5 m. The field is 
independent with respect to radius and depends on the z coordinate as B y (z) — 1.168 — 1.651 cos(27rz/0.5236) 
where the field is in Tesla and the distance is in meters. 



test 6 sector anticyclotron 
0.523598775 111111111 

1 

n s da q 

1.167764501 


10 -1.651468395 




1111 

0000000000000 
000000000000 



D.4 output for009.dat file 

The following shows the kinematic information of the reference muon and the three muons through one 
wedge sregion and one vacuum sregion. Each data line wraps around to form three actual lines. 

#test ring 180 MeV (H5b) 

# units = [s] [m] [GeV/c] [T] [V/m] 

# evt par typ fig reg time x y z 

Px Py Pz Bx By Bz wt Ex Ey Ez arc 

polX polY polZ 

2 1 0.0E+00 0.0E+00 0.0E+00 0.0E+00 
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O.OE+00 O.OE+00 1.8E-01 O.OE+00 O.OE+00 O.OE+00 1 . 0E+00 


10 2 1 O.OE+00 O.OE+00 O.OE+00 O.OE+00 

0.0E+00 O.OE+00 1.8E-01 O.OE+00 O.OE+00 0.0E+00 1 O.OE+00 


2 2 1 O.OE+00 O.OE+00 0.0E+00 0.0E+00 

0.0E+00 O.OE+00 1.8E-01 O.OE+00 0.0E+00 0.0E+00 1 0.0E+00 


3 2 1 O.OE+00 O.OE+00 0.0E+00 0.0E+00 

0.0E+00 O.OE+00 1.8E-01 0.0E+00 0.0E+00 0.0E+00 1 0.0E+00 


2 5 9.1E-11 O.OE+00 0.0E+00 2.4E-02 

0.0E+00 O.OE+00 1.8E-01 0.0E+00 -4.2E-01 0.0E+00 1 2.4E-02 


10 2 5 9.2E-11 7.8E-04 1.2E-04 2.4E-02 

1.1E-02 1.9E-03 1.7E-01 2.4E-10 -4.2E-01 6.4E-04 1 2.4E-02 


2 2 5 9.2E-11 1 . 1E-03 -7.4E-05 2.4E-02 

1.5E-02 -1.7E-03 1.7E-01 -6.2E-11 -4.2E-01 -4.1E-04 1 2.4E-02 


3 2 5 9.2E-11 7.8E-04 -1.6E-04 2.4E-02 

1.2E-02 -1.5E-03 1.8E-01 -6.2E-10 -4.2E-01 -8.8E-04 1 2.4E-02 




The orbit in the LAB (x, y) plane can be viewed with Historoot which can be obtained by downloading 
G4Beamline from www.muonsinc.com. With for009.dat as the input file to Historoot, the orbit in the 
LAB (x, y) frame is obtained by plotting (x + 0.5 m) cos(z/0.5 m) vs. (x + 0.5 m) sin(z/0.5 m). The radius 
is the LAB frame of the anticyclotron is (x + 0.5 m) where x is the the horizontal coordinate with respect 
to the ICOOL reference path which has a bent solenoid radius of curvature 0.5 m. Fig. [3] shows the orbits 
of three muons in this example ICOOL anticyclotron simulation which are stopped in about 200 ns. 
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ICOOL Example Muon Anticyclotron Orbits 
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Figure 3: Orbits of three color coded muons in an ICOOL simulated anticyclotron. The output is generated 
by Historoot. Each point in the orbit is at the location of the end of an sregion. Unlike the G4Beamline 
simulation, the magnetic field in this ICOOL simulation is independent of radius and does not include an 
inner magnetic bottle. 
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